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Abstract 



We present a coherent account of how the entanglement interpretation, 
thermofield dynamical description and the brick wall formulations (with 

the ground state correctly identified) fit into a connected and self-consistent 
explanation of what Bekenstein- Hawking entropy is, and where it is located. 
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1 Introduction 



The Bekenstein-Hawking entropy S'bh has been derived in so many ways 
and interpreted from so many points of view. All existing derivations agree 
that the Bekenstein-Hawking entropy is proportional to surface area A, 

It is an entirely superficial result, in the literal sense that it refers solely 
to surface properties. So either, a profound holographic principle is at 
work (all information about the black hole interior is somehow imprints 
on the horizon) [1 or existing derivations are "superficial", in the sense 
that they refer to an effective black shell entropy and don't probe the real 
black hole interior at all. Following this last idea, in this paper we present 
a basic model of black shell from the existence of thermal energy strongly 
concentrated near the horizon (thermal energy "wall") with respect to an 
uniformly accelerated observer, according to the equivalence principle. 

Among the different derivations of S'bh that have been proposed to pro- 
vide a microscopic explanation of the Bekenstein-Hawking Entropy [2l|3l|4], 
maybe the most promising and appropriate approach is that Sbh is entan- 
glement entropy, associated with observable and non-observable vacuum 
fluctuations correlated at the horizon with respect to an external observer. 
A program of this type was first clearly formulated by Bombelli et al. [S]. 
It was independently re-initiated by M. Srednicki [B] and by V. P. Frolov 
and I. D. Novikov [7]. From this approach the entanglement entropy pro- 
portional to the area of the dividing wall can be obtained. These striking 
results of entanglement entropy are entirely general, by no means confined 
to systems in thermal equilibrium. Accordingly, they do not within them- 
selves bear any clue as to the origins of the peculiarly thermal character of 
the Bekenstein-Hawking entropy. A comprehensive understanding of the 
Bekenstein-Hawking entropy requires a consistent blend of the entangle- 
ment interpretation with the thermofield-dynamical description [51 [SI [IH] • 
Moreover, this allows one to say not only what this entropy is but where 
it is located. The resulting picture (which we elaborate in Sec. 3) is of a 
"thermal atmosphere" extending a few Planck lengths above the horizon. 

On the other hand, the entanglement interpretation is implicit in, and 
is certainly closely related to the brick wall model introduced by 't Hooft 
[TTj . which will be shown below. In this sense, we use the modified brick 
wall model [H] and consider it as a model of black shell (a massive spher- 
ical shell compressed into a thin layer near its gravitational radius). This 
correct description (Sec. 3) shows that the thermal atmosphere is properly 
understood as excitations above an energetically depressed ground state 
(the Boulware state), and provides an accurate localization of these exci- 
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tations in the Hartle-Hawking state. Because this locahzation involves dif- 
ferential expectation values, the calculations are actually simplified, since 
the ultraviolet divergences cancel out. 

It is thus our aim in the following sections to provide a coherent account 
of how the entanglement concept thermofield dynamical description and the 
brick wall formulations fit into a connected and self-consistent explanation 
of what Bekenstein-Hawking entropy is, in terms of a model of a massive 
spherical shell compressed into a thin layer near its gravitational radius, 
considering the Hartle-Hawking state as an effective Hartle-Hawking state 
for a modified brick wall model 'TT . 



2 Thermally entangled Minkowski vacuum state 

Consider a real scalar field $ defined on the geometrical background with 
static metric given by 

ds' = ^f{r)dT' + ^+dz'+dy\ (2) 

for /(r) — r. 

One set of mode solutions of the corresponding wave equation for R : x > \t\ 
sector of Minkowski spacetime is the following: 

Mr, = Mr) e^^'^^+'^^y^ ] , (3) 

with a; = = (r, ,17 = ujk-^k^, where the modes ipn{r) satisfies the 
equation 

^Mr?j + [^^ - {k^ + kl + m") /(r)] ^r) = . (4) 

Now we introduce Killing-Boulware (KB)-modes Lp"^{x) for R and L : 
a; < |t| sectors 

' (x) = "^o (x, t) (x) , (5) 
where is defined in terms of the null coordinates u, w, 

Q,{x)^\{Q{-eu) + Q{ev)}. (6) 

According to thermofield dynamics, for a quantum description consider 
both quantization schemes, Boulware-Rindler and Minkowski, that is 

^{x) = X^n\^) = E ■ (7) 

e,f2 t,n 
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where, 



a[,^) = b^^^ cosh X - sinhx = e^^ fc^,^) e^^ , 



J>0) 



o 

(w>0) 

= e(iLj) e(5ee' Snn' ; 



with the vacuum states defined by 

«n'|0)M-0, bl^^\0)^^0, (^e>0). 
According to ([8]), these vacuum states are formally linked by 

From (fT3|) and ([T0| we can write 



(8) 

(9) 

(10) 

(11) 
(12) 

(13) 

(14) 



where, Izi)^^ 1 1^)31) excitations of the parochial Boulware states 
I and I , respectively, and 



n = {no , for all 57 with lu > 0} ; En — rifi ui , Z = 1 



n 

uj>0 



Also, from (IT4t we can define Boulware state \0)^^ over the complete space 
as I 0)^j depopulated of all Boulware modes \ n)gj^ in _R-sector, 

I (over i + i?) = Z-i E 1 2i) J )h • (15) 



2.1 Hamiltonian Formulation 

The physical sense of the transformation ^ is based in the invariant action 
and invariant Hamiltonian H under this transformation, according to 
the action 



sm = pjt (^E^^*^'(*)) > 



(16) 
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where the Lagrangian L is given by 

L = ^eL(^)($); ^ j &\^) £x, (17) 

{-ff°°'i>',t-(5'^''$,a<i>,6W$2)}e,(x), (18) 

with g EE detg^^. 

Then, from pT]) . we can obtain the on-sheh calculation for the expectation 
value of the local Hamiltonian operator in the global quantum state 

|o)m 

m(0 I 0)m = -Z-^ e-^^^S„ = (^)/3, (19) 

where the local Hamiltonian is given by 

Hi . / d^lL-W l-l i^n'' b^n^' + b'^P ■ (20) 
n 

The resulting ([T9|) means that a static observer in the region R perceives 
the vacuum state |0)^ as a thermally excited respecting vacuum state 
|0)jj. We can interpret it as the thermal feature of Minkowski vacuum 
for a local observer, due to his restriction to region R and the presence 
of an event horizon; whose effects are a shift of his ground state bellow 
the global Minkowski ground state and it perceived as thermally excited, 
respectively. That is, Minkowski vacuum is perceived by this restricted 
observer as excited thermally above his ground state. 

2.2 Minkowski thermal energy wall 

In order to show the existence of a thermal energy "wall" associated to 
the properties of the physical vacuum in the vicinity of the horizons, let us 
calculate the expectation value for the component Too{x,x') of the stress- 
energy tensor with respect to Minkowski and Rindler vacuum states for 
a scalar field. In general, for a scalar field, the expectation value of the 
stress-energy tensor Tai3{x) is given by 

{T^p{x,x')) ^V^p,W{x,x'), (21) 

where W{x, x') is the Wightman function and 

V^p. = a(„ dp>)-^ go^p {d-'dy + m^) . (22) 

Then, 

lim ( r»(x, x') )^_^ = d° d,. VFm_h- (23) 

T i-T' 
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where Wm-b. — Wm — W^r is the difference between Wightman functions for 
Minkowski and Rindler vacuum states, respectively. 

do do' Wm-r(x,x')|^,=^ = J2 effi^i „ i ^" ^o' 



UJ 



_ 1 



i'y'f2(£)r 



where the functions ififi{x) are 



then, 



n 

oo 1 



;/3|"l - 1 



Under WKB approximation, with 

the equation ([4]) leads to the differential equation 



with 



Then, we obtain 



do do> Wm-r = ^ 



where we approximate, for large k^^k^, 



l\k\ 







J2Y.O^ [Odp, also ^() = 2() 

fc, k ri=± 



(25) 

(26) 

(27) 
(28) 

(29) 
(30) 

(31) 
(32) 
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Moreover, k^^^{u!,r) was defined according to 

fc2(r;a;,fc„,J = 0, (33) 



tliat is 











pdp 




m 


/ 



kL.. + kL.. ^^-m^ + r'^p'. (34) 

(35) 





Thus, we finally obtain 

Jo eT - 1 {2Try 

where the extra term in — can be shown to be zero to the same approx- 
imation. 

This is the expected thermal expression for energy density of a hot scalar 
field, strongly concentrated near the horizon under WKB approximation, 
which is good near the horizon. 



3 Entanglement entropy model of black shells 

From the existence of thermal energy strongly concentrated near the hori- 
zon (thermal energy "wall") with respect to a uniformly accelerated ob- 
server, according to the equivalence principle, we can expect that identical 
thermal effects will occur, not only near the horizon of a very massive black 
hole but also near the exterior of a starlike object with a reflecting surface, 
compressed to nearly (but not quite) its gravitational radius. In particu- 
lar, the thermal description, summarized by expression (1361) . is the same 
as the one established for the modified brick wall model 12^. Just as the 
Minkowski vacuum, according to the description above, is explainable to a 
uniformly accelerated observer as a thermal excitation above his negative- 
energy (Rindler) ground state, close to brick wall a delicate cancellation 
between a large thermal energy and an equally large and negative ground 
state energy is manifested. Following this model we may approximate the 
total stress-energy (ground state -I- thermal excitations) (Tap)^, near the 
wall, to the Hartle-Hawking stress-energy (Tap)^^: 

(T^p)^^ « (T^p)^ - (T„/3)3 + (AT„^)„_ (37) 

i.e., effectively the Hartle-Hawking stress-energy, which is bounded and 
small for large masses. So, in very good approximation for this case, near 
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the wall we can use Hartle-Hawking state as an effective Hartle-Hawking 
state. 

In this context we use Hartle-Hawking state below, for the model of a 
massive spherical shell compressed into a thin layer near its gravitational 
radius. 



3.1 Thermofield dynamics of black shells 

We consider the generic situation of a real scalar field $ propagating on 
the geometrical background with static metric given by 

ds^ = g^g dt^ + gab dx" dx^ . (38) 

Such a field $ satisfies the field equation 

(□-m^)$ = 0, (39) 

where d'Alembertian operator □ is expressed by 

□ = {~9Y^^-do {{-9Y'^f'd.\ + {-gY^^-da {{-9Y'^g''''db\ , (40) 

with g = det 5^11,. 

In particular, for R sector in the Schwarzschild spacetime maximally 
extended {R : Z > \T\) with static metric 

ds' = -f{r)dt' + ^+r'dn^ (41) 

one set of mode solutions of Eq. (|39p is the following: 

^o(i,a:") = Mr) YirniO, ip) g-'"*, (42) 

where 

x = x"' — {r,d,ip) ,n = Lolm. 
Substituting Eq. (j42|) into Eq. ([39]) it is obtained 

There are two special solutions ^pn±{r) for the equation defined 
by the boundary condition 

^o±(r)«e±''^'-*, (44) 
for r — > ro (event horizon) 
or — )■ — 00 with defined by 
dv 

dr*^j^y (45) 
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which yield modes that are outgoing and ingoing modes on the horizon in 
R sector. So, = Loklm for A; = ±1. 

Now, in order to concisely write modes for R and L sectors {R : Z > \T\ 
and L : Z < |T|), consider t-modes ipn{x,t) for R sector and Killing- 
Boulware (KB)-modes ipQ\x) for R and L sectors: 

^n{x,t) = ^n{x)^=e-'^\ (46) 

v[^\x)^Mx,t)Qe(x), (47) 

where 

e,{x)^^{ei-eU) + e{eV)} , (48) 

Q is the unit step function, t runs backwards in L sector (e = —1) and 
U, V, are the Kruskal times (see Appendix A). 



With the K-G products 

i^Pn, ipn') = <^{^) Snn' , (49) 
(^ip^^\v^^?)^ee{uj)Snn'S,,>, (50) 

and 

e{uj) = sign(tj), 

Snn' = (5(a; - w') Skk' Sw Smm' , 
k = ±1. 

These (KB)-modes share the same static metric (HIT) but with different sign 
for the time coordinate. More precisely, (KB)-modes are positive frequency 
modes in Killing time Epsilon t in each of the sectors L, R. 

For a complete set of modes given by Eq. (j47| , we have the completeness 
relation 

I J2 M3L) vhW) =S'i3l- x!) , (51) 
n 

and the orthogonality relation 

J d^xj{x)(f*^,{x)(pn{x) = 6nn' (52) 
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where, 



, Imk. 



On the other hand, consider Kruskal-Hartle- Hawking (KH^)-modes Xn^ (x) 
which are positive frequency modes in Kruskal time U, V over complete 
manifold; i.e. Xn\^) is positive frequency in kruskal time if we > 0, and 
Xn\x) is negative frequency in Kruskal time if cj e < (See Appendix A). 
These modes satisfy the orthogonality relation 

(xL'\x[i<^) ="Hfeo'<5..', (53) 

and are connected with the modes ip^^^ (x) by the Bogolubov relation 

(54) 

where x is defined in terms of 

tanhx = e-'^l"l/''o, (55) 



XnHx) = coshx + (fin sinhx, 



and K„ is the surface gravity. From equations (|T7)) . ([M)) and (|170p it is also 
calculated 



/sinhx coshx 
2M 



(56) 



For a quantum description consider both quantization schemes, for 
Boulware and Hawking-Hartle, that is, for the modes (|Tf)) and ([55]) . re- 
spectively, 

Hx) = Yl xl^^i^) = E ^^\^) , (57) 
e,n e.n 

where, from Eq. ([54)) 



,M - hi-) 



b^' coshx — sinhx — e 



(58) 



and 



Wt^(-e) 



G = Gt^i^^e(w)ex6[,^'^&; 



^(+)t^(-) _^{+)5(-)t 



(w>0) 



(t.J>0) 



(59) 



Now, we postulate commutation relations 



= e(iLj) eJee' fca' ; 



(60) 
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with the vacuum states defined by 



0, b^^^\0)^ = 0, (0.6 >0). 



(61) 



According to Eq. ([55]). these vacuum states are formally linked by 



10) =e-*'^ 10! 



(62) 



This last relation needs to be treated with some caution, because strictly 
I 0)jj and 10)3 are unitarily inequivalent [9]. 

A more explicit form of the Eq. can be obtained by factorizing the 
operator exp{— iG} into its creation and annihilation operators, which is 
done by a generalization of the Baker-Campbell-Hausdorff (BCH) identity 
(see Appendix B): 



-iGn 



0] 



nn=0 



-5 ^ riQ w 



(+) (-) 



(63) 



where 



is the Boulware state with an equal number uq of 



correlated Boulware modes in the state of the field 51 ( i.e., it specified for 
the quantum numbers fl = uj,k^l^m with > 0) in the L and R sectors. 
Moreover: /3 was introduced in terms of Hawking temperature T^, such 
that 



Eq. (1551) was rewritten as 



(64) 



tanh X = e "0 = e 2 p"^ , 



uj>0. 



(65) 



and we have used 



cosh^ X = 



1 - e-P'^ 



-n p LU 



= Zn. 



(66) 



According to equations ((5^)) and the expression ([S^ can be reduced 
to 



10) 



n 

ui>0 



(67) 
(68) 
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where 



n = {no , for all O with io > 0}, 

En^Y.''^^^ ^ = ^e-^^^= n Za. (69) 

n n O 

ui>0 cj>0 

Then, the ground state on the full Kruskal manifold is a thermal entan- 
gled state. Vacuum state | 0)^ has been written as a global state, but it 
does not really exist as a global state defined over the L-\-R region, because 
(KB)-modes are not positive frequency in any globally defined regular time 
parameter. The state \ is empty of (KB)-modes (p^^ positive frequency 
in (future directed) Killing time t in each of the sectors L, R. 

Consider the parochial Boulware state \ defined over the complete 
Kruskal space: 

from expression it is clear that 

|0)H=^"^Ee"^'''"l^i)Bj2i)BK, (70) 

n 

where, Izi)^^ ^^'^ \lk)^Yi ^'"^ excitations of the parochial Boulware states 
I and I , respectively. 

On the other hand, From Eq. ([70|) we can think of | as | )jj depop- 
ulated of all Boulware modes \ n)-^^ in R sector. The parochial Boulware 
state I )3j^ is empty of positive frequency Killing modes ip^^'^ in the sector 
R: 

&[2^M0)3, = 0, OO. (71) 
3.2 Thermal energy of black shells 

Consider the expectation value for the component of the stress-energy 
tensor with respect to Boulware and Hartle-Hawking states for the scalar 
field defined over metric (pS)) . In general, for the scalar field $, the expec- 
tation value of the stress-energy tensor Tq,^ {x) is given by 

{T^p{x,x')) = V^p,W{x,x'), (72) 
{Tc,p{x)) ^ lim (TQ^(a;,a;') ), 

where 

To.fi = ^,c.^.f}-\ g^p $,7 +m2$2) , (73) 
for the minimal coupling; 

Vo^p, = dp')-^ go.p {d^dy + m^) , (74) 
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and W{x,x') is the Wightman function 

W{x,x') ^ {0 \^{x)<i>{x')\ 0) . (75) 

(j72p and (|75p become infinite wlien a;' — ^ x, because mode sum 
dcj (/?^(a;') iy3n(a;)|2:'=2: diverges as w oo. A regularization is required, 
and involves subtracting the (state-independent) Hadamard function 
H{x,x') from W{x,x') jl3]. But this wih not be needed here, because we 
are only interested in differences between Wightman functions associated to 
Boulware and Hartle-Hawking states, and H{x, x') cancels out of Eq. ([55)1 
and the corresponding {Tap). 

Regarding Boulware state | 0)^, the Wightman function is given by 

W^{x,x') = b(0 |$(a:)$(a:')| 0)b 

= b(0| ^^[.^VL^^(^)4'^VL'^*(^')|0)b, (76) 

where expansion (j57p was used. Thus, by resorting to the result 

B ( I b[i}^ I ) B = e(e L.) 5,,, Snn' , (77) 

we find 

W^ix, x') = ^n^(^) ■ (78) 

Since ip^^ {x) is proportional to 6e(x), given by Eq. (^5]). W^{x,x') — 
if X, x' are in opposite sectors L, K. 

Similarly it is obtained for | ) , 

W^(x,J) = h(0 |$(a:)$(a:')| 0)h 

^Y.^{,uj)^^\x)xt^\^). (79) 

By using the Bogolubov transformation (IMl) . W^H(a;, a;') can be expanded 



as 

{x, x') = ^ e(e c.) [cosh^ X vl^^ ^^n^* (^') + sinh^ x ^[T^' (^) ¥^[7'^* (^')' 



(80) 

for a:, x' in same sector, since ipQ\x) ip'q '^\x') oc Qe{x) 8_e(x') = 0, when 
X, x' belong to same sector. 
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From the expressions (|80l) and (1751) 
for a;, x' in same sector and where one used the identity 



(81) 



(W^-W^){x,x') 

1 



^1-1 



'^[p{x) ^[+)*(x') + ^l;\x) ^[;^*{x')] , (82) 



where it was introduced 



sinh X = 



tanh X 
1 — tanh^ X ' 

for tanhx defined by Eq. ([55t . 

Due to properties of ipQ\x) and (/?[^'^^(a;) it is necessary to restrict the 
calculations to sector R or sector L. Then, if we suppose x, x' belong to 
same sector, say R, finally Eq. ([82]) becomes 



(83) 



(Wh- VKb)(x,x') 

because in sector i?, 0_ = and 0+ = 1, 
then 

cp[-)(a;)^[->(a:')cxe_(a;)e_(x') =0 , 
(^[,+)(a:)^[+)*(x')(xe+(x)e+(x') =1. 
To calculate Tpp(a;, a;'), first we evaluate 

d„ 9„,(W^H - W^)ix,xX,^^ EE 5„, W^„_b(x,x')L,=, , 

a„a„,w^H-B(a:,a:')L,..- E e^J_i 

xa„a„. ^[+)(:r)^W*(a:') 



(84) 



(85) 
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where Eq. ([^^ was considered. 

The expression (|85|) can be reduced to 



duj ■ 



An Jq eP'^ - 1 
Y,{2l + l)\^^i{r)t 



by using 



Ifn (x) = ; (r) Yirn{d,(f) , 
m— — l 



(86) 



(87) 



47r 



= 2 / 

-oo Jo 



dio{). 



(89) 



The result (f86|) does not diverge since due to the difference between Boul- 
ware and Hartle-Hawking states introduced, the Bogolubov transformation 
([M)) supphes a convergence factor (e^''^' — 1)~^. 

To complete the calculation of Eq. ((86)) we need to known explicitly the 
function ipi^i. Thus, in terms of the metric (|4ip and introducing the new 
function ip{r) in the mode solutions (|87p such that 



f 



■ip{r) , 



the equation Eq. ([5^ leads to the differential equation 



rf2 



where 



1 (uj^ l{l + 1) 
7 IT~^^ 



ir'f)" 
2r2 



(90) 



(91) 



(92) 



and the event horizon is characterized by r = r^, /(r^) =0. 

In order to solve the equation Eq. (|9T|) we resort to WKB approxima- 
tion, then 



2n 



k{r) 



i k{r ) dr 



(93) 



where the new index 77 = signfc(r) corresponds to outgoing and ingoing 
waves: ery = ±1, which means that we need include 77 in collective index 
0, = ujlmrj. 
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Before continuing the calculation of the expression , it is important 
to consider the normalization of WKB mode ((93|) . So, spatial modes fn{x) 
are required to satisfy 

d^x'y{x) (pn'{x) (pn{x) = S^n' + Sqq,, (94) 

where 

f2 = ujlmri, Cl = — cj, Imrj, 
l{x) ^ V^{-g"") = .r^r^ sine. (95) 
In these terms, now we have 

Vuiriir) =— = i^uir^ir) , (96) 

de dip sin e ie,ip)Yim{e,ip) ^ Sw s^m' • (97) 

Hence Eq. (p4)) reduces to 

drf-^ru'i'n' (r) = {S{lo - oj') + S{lu + u:')} 6,„, . (98) 

Returning to the calculation of (15^ . the expressions ([M|) and (IMl) are 
substituted into Eq. to obtain 



where we approximate, for large I, 

^()« /Od/, also 5]() = 2(). 

Moreover, /„,^^(cj, r) was defined according to 

fc2(r;w,C,J = 0, (100) 

that is , 

r-' LUL.. + l) = j-m'- ^ ■ (101) 

'"^^ ^^"iwr ^ 2r^V7^Vu.(u. + i). (102) 
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Finally, substituting Eq. into Eq. ([M]) and multiplying by (-5™), 

we obtained 

-a'a,.w._..J^fl^^f. (103) 

Jo e- - 1 [2ny 

where the local proper energy per mode E has been defined by £^ = 
Then, /3w = 

Tir) = Tn^ , (105) 

which is Tolman's law, and 

Lo duj 

pdp= — p =EdE. 

Regarding the extra term in —T° arising from d'^ d^i + m? in Eq. (|74p . can 
be shown to be zero to same approximation. Then, 

lim(T;(x,x'))H-B =5«5„,1^H-B 

_ E A-Kp^dp 

(106) 

This is the thermal expression for energy density of a hot scalar field cor- 
responding to eq. (p6|) . 



3.3 Entanglement entropy of black shells 

To calculate the entanglement entropy 5* consider the density matrix 
obtained from Eq. (j70p and the reduced density matrix p'--^\ given by 

Ph = |0)hh(0|, (107) 
= trs,-, Ph , (108) 

where the trace is taken over the degrees of freedom corresponding to sector 
L. 

Thus, we may calculate entropy as 

S' = -tr(p^+) lnp(+)), (109) 



16 



with the reduced density matrix 

p(+)=Z-i^e-'3^- n^+A^Jn^+^ , (110) 



where Z is expressed by (j69D , and 

OO - 

E ^""'^ = • (111) 

The thermal system described above and associated to entanglement 
entropy actually allow us to think that entropy arises physically located 
near the horizon. Thus, we can calculate it from the expressions (|69p and 
(|llip in terms of the partition function Z: 

S = -(3—lnZ + \iiZ , (112) 

which corresponds to 

S = l3U + liiZ, (113) 

with 

u = z-'Yl ^li^"'''" = ^ • (^^^) 

n 

Since, for any function f{uj) which goes to zero as a; — )■ oo, one can 
write, 

V /(..)= / dujN{u:)f{uj), (115) 

"'0 



then 



POO 

lnZ = J2 In Zn = ^ /{lo) = / dtu N{lj) /(w) , 

o o "'0 



(116) 



where 



/(-) = In [r^-p^ ) ' (117) 

and N(io) duj is the number of modes ^n{x°') which fall in the range (cj, w + 
duj) for all admissible fc, m. 
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In other words, in order to calculate S associated to given by Eq. (jlOSp . 
all what we need is to find N(uj) finite, i.e., count modes: 

lnZ= Y: (2^ + 1)1i^(y3^) ■ (118) 

Hence, to get S finite we need to restrict the discrete sets = 0, 1, 2, • • • } 
and {n = 1, 2, • • • }, which for the shell model is satisfied by the follow- 
ing restriction to the sets {?} and \n\'. the discrete set \r\ ~ ±1; Z = 
0, 1, 2, • • • ; m = — Z, • • • , is restricted by the condition for the radial 
wave number fc, defined by Eq. (|92|) . such that 

fc2(r;cj,/) > 0, (119) 

that is, k is real; otherwise '^(r) decays exponentially and is effectively 
zero. Hence, for a given w, N(uj) is finite or can be made finite by suitable 
boundary conditions. 

In that sense, consider a finite-size confinement for the entanglement sys- 
tem, according to brick wall model, with Dirichlet boundary condition 
ipif) = at inner boundary r — + e. 

WKB solution of the differential equation (I^TI) is the following expres- 
sion 



V'i„(r) = sin / fc(r'; toin, I) dr' , (120) 

y'k(r;UJin,l) Jro+e 

where for an admissible mode, both k{r) and ipir) must vanish simultane- 
ously at the n*'* node r = rn- 

k{r';uJin,l)dr' ^nn, (121) 

'ro+e 

k^{rin,uji„,l) = 0, (122) 
k\r';uji„,l)>0,r' <rin, (123) 

which means that fc(r) is large near + e and decreases with increasing 
r. If > [j2^ then fc^(r) must eventually become zero; say, after n 

oscillations of ■0(f'), r ~ r„. 

2 — — — 

Since in Eq. ([M)) is very large near the inner wall, high-Z allowed modes 
are very numerous. So these modes are confined entirely to a thin layer 
near the inner wall, with the contribution 

(124) 
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This dominant inner-mode contribution to In Z and S is proportional to 
r^. On the other hand, outer modes (which satisfy boundary conditions at 
outer wall r = R) give a contribution proportional to R^. This basically 
explains why the detailed calculations below gives 

S = 5„.„ + , (125) 

where 

S^^u (X horizon area. (126) 
From Eq. (jllSp we can write 

lnZ = ^(2/ + l)ln(^— i^) , (127) 

l,n \ ' / 

where it was considered that Dirichlct boundary condition is imposed at 
inner boundary, o alternatively Neumann boundary condition, which does 
not affect counting. 

Since I, n are integers, the sum over / and n can be replaced to a good 
approximation by a double integral. Instead of ^, n as independent variables 
of integration, we switch to l,uj as independent variables with n = n{uj,l), 
and use the expressions (|121|) and (|122p to write 



1 



n{uj,l) = - I k{r';uj,l)dr' ;k'^{r{uj,l);uj,l) =0. (128) 

Jro+e 



dn{iu,l) 1 r("^') dk{r'-LO,l) , , 

— ' ■ dr 



duj n Jro+e du) 



l^lpllkir{u:,l);u:,l), (129) 



where the second term is zero by Eq. (I128p . 
Now, from the equations (IT?7)) and 

\nZ^ - III dldujdr' (21 + 1) 



TT 

(fe2(r';L;.'/)>0) 

where the limits of integration are the surfaces k'^{r'; iv, I) — 0, r' — + e, 
r' — R, with / > and u > restricted only by the condition k'^{r'; to, I) > 
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0. 

Integrating Eq. (|130p by parts, it is obtained 
1 



InZ : 



dr' {2l + l)dldujk{r';uj,l) 



/3 



- 1 



(131) 



On the other hand, 

J {2l + l)dlk{r';uj,l) 

where 



1 



Lmax(tJ,I-') 



1 2 



'773 



(132) 



for given by the vahie that makes fc^ = 0; i.e., from Eq. 



2r'2 



f{r') 



Substituting Eq. (fT32)) into Eq. (IT3T|) 



In Z = — I duj 

TT 



1 



_ 1 r'T/ 3 



(133) 



(134) 



Prom Eq. (I134p we can obtain usual thermodynamic expressions in 
terms of the Helmhohz free energy F and the average energy U: 



where 



N{uj) 



37r 



U 



F ■ 
^ r'^dr 



InZ 



N{uj) duj 
e^" - 1 ' 



Pir) 



Lu — \ m + 



2r'2 



fir) 



(135) 



(136) 



dp 

1 



dujN{uj) 



NiLj)~-^{ujN{u)) 



-R -ujN'iu). 



(137) 
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with N'{uj) defined by 



Then, entropy S gives 

s^m-n^ff"-^^^l±^^.. (139, 

To convert the expressions above to a statistical-thermodynamic form, 
can change variables of integration from r and a; to r and p, where 



we can 

2 



P-^-™-^. ("0) 



fir) 2r2 

Then, 



pdpdr — J, dr . (141) 



In these terms, the expressions (|137p and (|138p become 

C,./" 4„.„. rf£41l£!iE, (142) 

where we have exphcitly restored h — 1 and introduced the locahy measured 
energy of mode with frequency uj, E — with uj(j), r) given by Eq. (|140p . 



and the effective local momentum p 
Then 



Thus, 



with 



nr) = ^- (143) 



U = in r'^ p{r) dr , (144) 

J ro+c 

f°° E Anp^dp , , 

= / , (145) 
Jo eT — 1 n. 

-B^ = p2 + + « + . (146) 



2r2 

Finally, from Eq. 
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Now, resorting to the equations (|136l) . (|140l) and (|14ip 



^ p2 f A 2 dr J. f p'^e'^'^ A-Kp^ dp 



Thus, 



where 



5= / 47r r2—s(r), (149) 
1 p p^ef 47rp^dp 

The analysis above corresponds to brick wall model [TTJ [T^ . According 
to this model, the integrals (|144p and (|149p are dominated by two con- 
tributions, for large r = R and for small + e. The former corresponds 
to a volume term, proportional to |7rr^, which represents the entropy and 
energy of a homogeneous quantum gas in a flat space at a uniform temper- 
ature . The latter is the contribution of gas near the inner wall r = R„. 
Then, for this last contribution is required to introduce the ultrarelativistic 
approximations 

47V 3 37V 4 
s = ^r^/9=— T. (151) 

Substituting Eq. (|15ip into Eq. (|149p . the wall contribution to the total 
entropy is obtained: 

iV 1 
90^4' 

where A'' accounts for helicities and the number of particle species, A is the 
wall area and a is the proper altitude of the inner wall above the horizon 
of the exterior geometry. 

Now, depending on a, we can obtain the Bekenstein-Hawking entropy 
from Eq. 

Swall — 5bh, (153) 
where a has been adjusted by invoking quantum gravity effects. 



Swall = -TTT. 5-7^, (152) 



4 Conclusion 

We have presented an integrated and detailed explanation of how the en- 
tanglement interpretation, the thermofield description and the brick wall 
formulas (properly interpreted as referring to thermal excitations above the 
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Boulware ground state) fit together to form a coherent, self-consistent ex- 
planation of what Bekenstein-Hawking entropy is, and where it is located. 
Nearly all the other brick wall papers are just formal calculations which, for 
instance, give no clue why the large entropies they derive near the horizon 
should sit on a region which is empty and nearly flat. 

Really, in this paper we have introduced a model of black shell in terms 
of an effective Hartle-Hawking state, which is externally indistinguishable 
from a Schwarzschild black hole with respect to its thermodynamic prop- 
erties. For this reason, the model considered above can be interpreted 
as an effective calculation of the entanglement entropy associated with a 
Schwarzschild black hole, i.e., S'bh can be considered as entanglement en- 
tropy, which is well defined near the horizon and presents a thermal nature 
according to thermofield dynamics of black holes. 

On the other hand, by using thermodynamical arguments we suggest 
that the brick wall model might be considered as a model of black shell 
in order to get an operational approach to black hole entropy [13]. The 
topped-up Boulware state (TUB) defined there may be called a general- 
ized Hartle-Hawking state more than an effective one. Indeed, it becomes 
the Hartle-Hawking state in the limit when the shell approaches its gravi- 
tational radius. 

The most interesting interpretation with respect to the model of black 
shell developed above is just that maybe all existing derivations of S'bh are 
"superficial" , in the sense that they refer to an effective black shell entropy 
and don't probe the real black hole interior at all. The viewpoint advo- 
cated in this paper suggests that it is possible to entertain the suspicion that 
all derivations of the Bekenstein-Hawking entropy formula, which differ so 
vastly in appearance, are just disguised variations of the same derivation. 
All derivations lead to the same formula because all calculate entropy of 
the same object: a black shell instead of a hole. 
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A Positive frequency modes 
A.l Definition 

Let ln_(_ X be defined, for a real number x, by 



ln+ a; = In |a;| + , —oo < x < oo , (154) 



where 

e{x) = sign (x) . 



Then, 



e 



(155) 



are positive frequency functions in a; ( for both signs, ±a, a real), i.e. 

A±(w)e"*""'dcj. (156) 



g±za lii+ X 



This assertion, whose proof is shown bellow, it will be taken as definition 
of a positive frequency function. 

A. 2 Extension to an analytic function 

To prove the assertion above, extend the expression (|154|) to an analytic 
fuction ln+ z by defining: 



ln_|_ z 

that is. 



real on lower imaginary axis 
branch cut in upper half-plane , 



ln+z = ln|2;|+i(argZ+-), — ^ < argz < ■ (157) 

Then ln_|_ z is regular in lower half-plane. 

g±iain+ z regular and bounded in lower half-plane. 

Hence 

A±{uj) = —I e±'"'"+^e^"^dz = if w < . (158) 



— oo 



Finally it is defined 



Then, 
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In^ a: = In + —e{x) e , — oo < a:: < oo ; e — ±1 . (159) 



±min, X J positive frequency in a; for e = +1 
6 ^ are s 

I negative frequency in x for e = — 1 . 
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A. 3 Application to maximally extended black hole ge- 
ometry 

Let sectorial functions be defined by 



ee{x)^-{ei-eu) + eiev)} , 

where O is the unit step function and U, V the Kruskal times. 
U and V can be defined by 



T±Z : 



dV 



dU 



\duj 



dt± 



dr 

W)' 



Consider the definition 

2k„ t+ = ln+V - ln+ ?7 = In 



ITT , 



+ y(e(V^)-e(f/)) 



In 



+ i7r(e+-e_), 



and the extension 



= In^ y - lue U = \vl 



+ -(e(F)-e([/))e 



In 



+ i7r(ee - e_e) , 



then 



e=^*"*' are • 



positive frequency functions in (7, V for e = +1 
negative frequency functions in (7, V for e = — 1 , 



(160) 



(161) 
(162) 



(163) 



(164) 



A. 4 Useful relations 



(1) e, + e_, = l, Oe - e_, = - e {e(y) - e(f/)} 



(2) e- 



= e 



I e 



Oe + e = "0 ' e_ 



(3) Let e' = e(a;), then 



1 ttI^I _i 7r|u>| 



(165) 
(166) 

(167) 
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Defining 

X = X{^) by tanlix = e "o , (168) 
Eq. (|167p can be written as 



sinli X / V cosli X 



(169) 



[(sinhx)(coshx)]'''" e"'""''-'' = =e 

X [coshx0c + sinhxB-e] , 



(170) 



B Baker-Campbell-Hausdorff identity 

If two operators A,B and their commutator C satisfy the commutation 
relations 

[A,B]=C, [C,A]^2n^A, [B,C]^2n^B, (171) 
for some number n, real or complex, then for any parameter x, 

An outline of the derivation of the identity (I172p can be shown by 
introducing the operator [15) 

F(x) = e^^^^)'^ e^(>^)^ e^^'^)^, (173) 

where x, y, z are undetermined functions. The idea of the proof consist in 
choosing these functions so that F is reducible to e^^^'^-^\ 
Differentiating Eq. (I173P with respect to X; 

F'^F'ix) = x'e-"^(e-f ^yle''^)e"^ + i/'g-^'^Be"'^ + z'C, (174) 

we obtain 

F-'^F'ix) = x'(e-2'"^A + yC - y^me^'^^'B) + y'e^^^B + z'C, (175) 

where 

g-aB^paB ^ ^ _,_ ^] [[^^ B], B] + ... 

= A + yC — y^mB, 
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considering 

[C, A] = 2mA , [B, C] = 2mB. 

Now, it is required Eq. (jl75p to equal A + B. Equating coefficients of A, B 
and C yields three equations for x,y,z : 

x'e'^"'' = l, {y' ~mx'y')e^"'' ^1, z' + x'y^O. (176) 

If we define /(x) = e^^™^, the first and third of these equations (|176l) then 
give x' = /^^, y = f /2m. Substituting into the second equation results in 
a second-order linear equation for /^/^. Finally, in order to find the general- 
ized identity (jl72p . consider the solution of this last linear equation, subject 
to the initial conditions x«?/«X7-2^~0 when % — )■ is f{x) = cosh^ nx- 

Thus we can define 

An = 0^ , Bn = -Oj^ , uj >u. 

Then, 

On -b^ 0^ + On 0^ , 

and the commutation relations (|17ip are satisfied with n — 1, which means 
that Eq. (|172p is applicable. 

According to Eq. ([5^ . it is directly obtained 

e"*^ \ 0)^ = eA^"+^"^ 10)3 

1 ,(tanhx)!)y"'+fcL"' \n\ 



coshx 

00 

— ^(tanh"x) n^n\n^n^) , (177) 



cosh V 



which corresponds to Eq. (|63 
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